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Abstract 

We have studied the InxGai-xAs/IuyAli^yAs (001) interface using first-principles ab-initio pseu- 
dopotential calculations, focusing on the effects of alloy composition and strain state on the elec- 
tronic properties. In particular we estimate a valence band offset (VBO) of 0.11 eV (InGaAs 
higher), including spin-orbit and self-energy corrections, for a strain-compensated configuration 
with homogenous composition x = y = 0.75 on a lattice-matched substrate. Unintentional compo- 
sition fluctuations which are typically limited to a few percent and different short-range order effects 
give rise only to small variations on the VBO, of the order of 0.1 eV or less, whereas intentional 
substantial changes in the alloys composition allow to achieve a high tunability of band offsets. 
We predict a VBO varying in a range of about 1.1 eV for interfaces between the pure arsenides 
in different strain states as extreme cases of composition variation at In2.Gai_a;As/InyAli_yAs 
heterostructures. 

PACS numbers: PACS: 73; 73.40.Kp; 73.20.-r; 68.35.-p 
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I. INTRODUCTION 



GaAs, AlAs, and InAs form the family of common-anion III-V conventional semicon- 
ductors covering the widest possible range of energy gaps, apart from nitrides.— They are 
therefore particularly suitable to be combined into alloys to form In3;Gai_a;As/InyAli„yAs 
hetero junctions whose electronic properties can be easily tailored according to the techno- 
logical needs, acting on composition to control and intentionally modify the valence and 
conduction band offsets (VBO and CBO), namely for band-offset engineering •^''^'^ 

The use of alloys in heteroj unctions has also some drawbacks. Beside controlled vari- 
ations in the average composition, unintentional composition inhomogeneities could be 
present in epitaxially grown alloys and heterostructures, as detected by tunneling micro- 
scope techniques. -^"^ Their origin can be ascribed to several mechanisms: inhomogeneous 
incorporation of the alloy components during growth, atomic diffusion at the surface during 
growth induced by strain inhomogeneities arising from stress relaxation and/or interface 
roughening, and also post-growth atomic interdiffusion with or without thermal annealing. 
Although in high-quality alloy-based nanostructures and devices such inhomogeneities are 
minimized, their residual occurrence can affect the electronic and optical properties . 

Another source of variations of the electronic properties in alloy-based heteroj unctions is 
the occurrence of spontaneous ordering in the constituting alloys. In the last years consid- 
erable theoretical and experimental effort have been devoted in investigating the effect on 
the band alignments and related properties . ^^i^'^i^^i^^ 

It is important for device design not only to predict the value of band offsets at hetero- 
junctions with given composition, but also to estimate the effects of both composition fluctu- 
ations and ordering. We address this problem here, focusing on the Ina;Gai_a;As/Inj,Ali_j,As 
heteroj unction and studying by accurate ab-initio simulations: (i) the band offsets at a 
given nominal composition x = y = 0.75, presently of particular technological interest,— 
considering different structural configurations; (ii) the effect of possible realistic composi- 
tion fluctuations with respect to the nominal one. 

A computational approach that has been used to simulate several alloy-based hetero- 
junctions, such as AlGaP/GaPi^, GaAlAs/GaAs^^ and AlGalnAs-based heterostructures^^, 
is the virtual crystal approximation (VGA). This approach considers an average crystalline 
field acting on the electrons rather than the effect of the individual atoms and therefore is 
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particularly valid where the elements forming the alloy are very similar, such as the Al and 
Ga atoms in the examples above. By its nature, the VGA has intrinsic limitations in the 
atomic-scale description of the materials (e.g. it cannot describe lattice distortions). Fur- 
thermore, it has to be applied with care also in the calculation of macroscopic properties. 
It is known for instance that it gives positive deviations from the experimentally observed 
linear behaviour of the alloy lattice parameter with composition (Vegard's law) which is 
instead recovered using a fully atomistic approach, with "real" instead of "virtual" atoms. 
To overcome this limitation, the VGA could be applied to calculate the average electronic 
properties only after having considered the correct structure.-^^ In a previous work we have 
used the VGA for a preliminar study of the electronic properties of In^Gai_^As/Inj^Ali_j^As 
heterojunctions.3. Here we study more accurately the case x = y = 0.75 using a fully atom- 
istic approach with "real" atoms. This is mandatory in order to investigate the effects of 
realistic composition fluctuations and ordering.— 

We predict a VBO of 0.11 eV (InGaAs higher) and we estimate that realistic compo- 
sition fluctuations and ordering effects are small and not exceeding s^O.l eV. For com- 
pleteness, we have also estimated the maximum range of tunability of the offsets at 
In^.Gai_i.As/Inj^Ali_j^As heteroj unctions by varying intentionally the composition of the 
alloys with variations of x and y separately up to the limit of pure binary semiconductors. 

In the next Section we present our computational approach. Section III is devoted to the 
preliminar study of the three binaries (GaAs, AlAs and InAs) constituting the junctions, 
both in their bulk phase and in strained configurations, and of the two Ino.75Gao.25 As and 
Ino.75Alo.25As alloys in strain-free configurations. In Sect. IV we report the results for 
the Ino.75Gao.25As/Ino.75Alo.25As heterojunction in the hypothesis of a uniform homogeneous 
composition and we discuss the effects of different ordering and composition fluctuations. 
We complete the work considering in Sect. V the interfaces between the binaries, in different 
possible strain configurations. Finally, we draw our conclusions. 

II. COMPUTATIONAL DETAILS 

We perform state-of-the-art first-principles non-relativistic pseudopotential self-consistent 
calculations within the Local Density Approximation (LDA) to Density Functional Theory 
(DFT) using the PWSGF code.-^ We address the reader to Appendix A for details con- 
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cerning the inclusion of spin-orbit corrections a posteriori. Periodically repeated supercells 
are used to simulate both the bulk alloys and the heterojunctions. We use norm-conserving 
pseudopotentials with d electrons of Ga and In in the core and the non-linear core correction 
for In. This choice is very convenient since-compared with a complete neglect of the effects 
of d electrons-with a very limited computational effort it improves significantly the descrip- 
tion of the structural and electronic parameters, as it is widely reported in the literature 
for the traditional III-V compounds. We address the reader to Appendix ^ for a deeper 
discussion based on additional accurate all-electron calculations on the combined effects of 
p — d coupling and spin-orbit and/or strain splitting in these systems. 

The plane wave basis set is expanded up to a kinetic energy cutoff Ecut = 20 Ry which 
gives well converged equilibrium structural and electronic properties for the systems consid- 
ered. Test with 25 Ry of cutoff have been done for the binary and ternary compounds. The 
integration over the Brillouin zone is performed using the (4,4,4) Monkhorst-Pack mesh for 
the FCC cell and corresponding meshes for the various supercells. 

For all the heterojunctions we consider pseudomorphic growth conditions, i.e. we fix 
the in-plane lattice-constant a\\ and the possible lattice mismatch between the constituting 
materials is accommodated without defects by a strain field e. Far from the interface the 
strain field is uniform and can be described by the perpendicular lattice parameter a±. The 
macroscopic theory of elasticity predicts:— 



where is the cubic lattice parameter of the strained material, cn and cu are its elastic con- 
stants, and e|| and e± are the components of the strain in direction parallel and perpendicular 
with respect to the substrate. The use of Eq. (1) with Cij calculated by first principles to 
determine a± is equivalent to a direct optimization of the atomic positions along the growth 
direction by stress and total energy minimization, but it is more convenient since a change 
of the substrate does not require a new self-consistent calculation. For the heterojunctions 
we take the prediction of Eq. (1) as a starting guess and then we perform a further opti- 
mization of the atomic position in order to obtain the correct interatomic distances also in 
the interface region, where they can differ.— We optimize the internal supercell structure 
until forces are less than 0.02 eV/A. The effects of these final structural optimizations on 
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VBO amount up to :^0.1 eV. 

As explained in details in Ref. 0, we calculate the VBO by summing two contributions: 
VBO = AE^ + Ay, i.e. the band structure term AE^ (the energy difference between the 
relevant top band-edges of the two materials in their appropriate strain state, measured with 
reference to the average electrostatic potential in the corresponding bulk crystal), and the 
electrostatic potential lineup AV containing all interface-specific effects and extracted from 
the self-consistent charge distribution in the heteroj unction supercells. We report a positive 
value for VBO (CBO) at the A/B interface if the valence band top edge (conduction band 
bottom edge) in A is higher than in B. Our final estimates of VBO(CBO) include also spin 
orbit effects added a posteriori using experimental data (see Appendix A). 

III. BULK CONSTITUENTS 

A. Binary semiconductor compounds 

Table H] summarizes the calculated relevant equilibrium structural and electronic prop- 
erties of the bulk binary compounds. The agreement between the experimental and the 
calculated lattice constants is within which is acceptable since the mismatch between 
InAs and GaAs (AlAs) is much larger (~7%, see Table ITT| . The theoretical elastic constants 
Cii, Ci2 and bulk moduli are equal to the experimental values within a few % for GaAs and 
AlAs, whereas the discrepancy is larger (by approximately 10-20 %) for InAs. Incidentally 
we note that the value of bulk moduli calculated directly via the MurnaghanSI equation of 

Cii + 2Ci2 

state from the total energy curve satisfies the well known relation: B= within the 

3 

numerical accuracy (estimated to be 10 Kbar). Therefore, our choice of pseudopotentials 
is justified and we can afford the study of the alloys and of their heterojunctions on a realistic 
ground, with the possibility of a correct description of internal distortions. 

The energy gaps are systematically smaller than the experimental values, as typically 
occurs in DFT-LDA, the worst case being AlAs. 

B. Ternary semiconductor alloys: Ino.75Gao,25As and Ino,75Alo.25As 

We study here the bulk properties of the pseudobinary semiconducting alloys which are 
constituting the heterostructure at the nominal composition x = y = 0.75. For definiteness 
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we focus first on the Ino.75Gao.25 As alloy. 

We consider only three different very simple ordered structures corresponding to the 
most homogeneous distributions of the different cation types (i.e. on the smallest possible 
scale compatible with the composition): Luzonite (labelled "L", with 8 atoms unit-cell), 
Famatinite ("F", with 16-atom unit-cell) and (001) 1+3 superlattice ("S", with 8-atom unit- 
cell). In the first one the Ga atoms are arranged in a SC lattice; therefore this structure 
has the same cubic symmetry which is typical of the pure binary semiconductors.— In the 
second structure the Ga atoms are arranged in a tetragonal body-centered structure, and 
their relative distance is greater than or at least equal to one cubic lattice parameter. The 
third structure is anisotropic, and it is the only one (among the three considered) where Ga 
atoms are next nearest neighbors. Larger supercells or alternative approaches accounting 
for compositional disorder on larger scales would be necessary for a complete treatment of 
the possible effects of randomness and for statistical analysis , but this goes beyond 
the purpose of the present work. 

Although the structures considered here are very simple, they allow us to catch the main 
structural features of the alloy, that we summarize in the following. 

The calculated equilibrium average lattice parameter is aaiioy=5.86 A for all the three 
structures, almost equal to the Vegard's value. 

In each structure, the nearest-neighbor (NN) pairs are of two types (Ga-As and In-As), 
with an occurrence proportional to the corresponding cationic concentration; the calculated 
individual NN distances (see Table IIIIll show a bimodal distribution centered around two 
values quite close to the bulk-like Ga-As (2.40 A) and In-As (2.58 A) values, as typically 
observed in most pseudobinary semiconductor alloys:-^- more precisely, the Ga-As distances 
in the alloys are within the range 2.41-2.45 A, whereas the In-As distances are within 2.56- 
2.60 A. The weighted average of all the anion-cation bond lengths is 2.54 A, coinciding with 
the NN distance in the ideal unrelaxed structure. 

The pattern of the next nearest-neighbor (NNN) environment is more complex, due to 
the presence of both As-As and cation-cation pairs. The anionic sublattice is rather dis- 
torted with respect to the ideal zincblende structure, due to different types of NN cations 
surrounding As, whereas the cationic sublattice is only slightly distorted, having only NN 
As atoms. As a consequence, the NNN distances of the cation-cation pairs (In-In, Ga-In 
and Ga-Ga)^ are quite close to the corresponding common value in the ideal undistorted 
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structure, i.e. 4.14 A: they are exactly equal to this value in the L and F structures, and are 
within a range of 0.15 A (from 4.07-4.22 A) in the S structure. The As-As NNN distances 
are rather scattered with respect to the ideal undistorted value, even in the L structure 
where the high symmetry limits the internal distortions: the As-As NNN distances have a 
quite broadened distribution in a range of about 0.44 A, from 3.92 A to 4.36 A. 

The Ino.75Alo.25As alloy shows structural properties very similar to the one with Ga, as 
it can be seen again from Table ITTTl Neglecting the tiny mismatch between GaAs and AlAs, 
it is perfectly lattice-matched with Ino.75Gao.25As. 

The electronic structures of the alloys in the three different structures are rather similar, 
as shown by the comparison of their band structures and of their density of states in Figs. 
and 121 with some differences. In all the three different structures considered, the calculated 
gaps are slightly smaller than the linear interpolation between those of the parent end- 
points, suggesting a small positive bowing in qualitative agreement with the experimental 
data available:— within LDA, we calculate Eg in the range 0.49-0.55 eV for Ino.75Gao.25As 
and in the range 0.82-0.90 eV for Ino.75Alo.25As (which has a direct minimum gap at this 
composition), to be compared with the linear interpolations 0.65 eV and 0.90 eV respectively. 

Although we do not perform an exhaustive study, we can however give a rough estimate 
of the effect of the structural order on the energy gap from the variations of the calculated 
values for the different ordered structures considered here: the maximum calculated variation 
of the gap is 0.08 eV, small but higher than our relative numerical accuracy which amounts 
to 20-30 meV. Incidentally, we notice that a comparable effect (0.06 eV) has been determined 
from photoluminescence excitation spectroscopy as band gap difference between the ordered 
and disordered epitaxial Ga^Ini_^.P.^^ Higher band-gap reductions (up to 0.25 eV and 0.18 
eV for Ino.5Gao.5As and Ino.5Alo.5As respectively) are instead predicted theoretically for a 
fully disordered— i'CuPt-[lll]-ordered transition, which is accompanied by the formation of 
important piezoelectric fields:^^ however, in real samples the degree of spontaneous order 
is not perfect, and the effect is expected to be <0.1 eV. We also notice that a band gap 



reduction up to 0.06 eV is predicted in Ref. 



12 in epitaxial Ino.5Gao.5As alloy, due to lateral 



composition modulation (which can be also seen as a ordering effect, but on a different 
length scale) with respect to an homogeneous configuration. 
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IV. THE /no.75Gao.25^s/-/'?^o.75^^o.25^s HETERO JUNCTION 



We consider the heterojunction between the two alloys Ino.75Gao.25As and Ino.75Alo.25As 
on a lattice-matched (001) substrate. This heterostructure indeed has been successfully 
grown in pseudomorphic and almost unstrained configuration by inserting InAlAs buffers 
with graded In composition on GaAs substrates . i7i35^36^37^38 ygj^y recently it has been reported 
on the achievement of a two-dimensional electron gas in these quantum well structures, with 
reduced carrier density and high mobility.— 

We recall here that also other particular compositions are of technological interest. 
In^Gai_^As/Inj^Ali_j,As hetero junctions grown with an In concentration x ~ 0.3 on an un- 
strained metamorphic buffer are widely used in microwave and optoelectronic devices.— 
Devices based on Ino.53Gao.47As and Ino.52Alo.48 As alloys are also developed for a wide variety 
of optoelectronic and high-speed electronic applications, since for these particular composi- 
tions they can be easily grown lattice-matched to the commercial InP substrate. 

We simulate the Ino.75Gao.25As/Ino.75Alo.25As interface by using tetragonal supercells 
made of three slab unit cells for each constituent alloy along the (001) direction. For the sake 
of definiteness we consider the same structure and configuration for the two alloys which 
form the interface. 

Since the alloys in the L and F structures are very similar in terms of structural and 
electronic properties, we consider for the calculation of the VBO only the L and the 5* 
structures. For the L structure the heterojunction supercell contains six simple cubic cells for 
a total length normal to the interface of 6ao=35.16 A, with oq the average lattice parameter of 
the alloys. The cation profile along (001) is: . . .X — In—X — In—Y — In — Y — In — . . ., where 
X indicates the mixed-cation planes with 50% Ga and 50% In, and Y refers to the mixed- 
cation planes with 50% Al and 50% In. We report in Fig. El upper panels, the macroscopic 
averaged of charge (solid lines) and potential (dashed lines). The calculated VBO with the 
optimized atomic positions and with the spin-orbit effects included a posteriori^*^ is 0.07 
eV, with a numerical error of the order of 20-30 meV. Self-energy corrections to the valence 
band top edges have to be considered for a more realistic estimate. From the values given 



in Ref. 



42, properly scaled to be adapted to our calculations,— the self-energy correction to 



our IDA VBO is 0.04 eV, thus giving a final VBO=0.11 eV. 

The VBO obtained from supercell calculations can be compared with the linear inter- 
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polation between the parent end-points: the GaAs/AlAs interface {x = y = 0) and the 
trivial case of the InAs/InAs homojunction [x = y = 1). The calculated LDA VBO at 
the unstrained GaAs/AlAs interface, including spin-orbit effects but not self-energy correc- 
tions, is VBO=0.44 eV, which well compares with previously reported values.^ The linear 
interpolation for x = y = 0.75 gives 0.11 eV, to be compared with the supercell calculation 
giving 0.07 eV. We thus confirm in this work — in sign and magnitude within the numerical 
accuracy — the deviation from linearity that was roughly estimated in Ref. Q with slightly 
different technical details and composition. 

In order to estimate the possible effects of short-range order/disorder, we then consider 
the case where both alloys are simulated with the S structure, which is a limiting case of 
maximum ordering in a particular direction, here considered as the growth direction. In this 
case the total length of the supercell normal to the interface is equal to 12ao = 70.32 A. The 
cation profile is the following: . . . — In — Ga — In — In — In — Ga — In — In — In — Al — In — . . .. 
The charge and potential profiles are shown in the lower part of Fig. IHl Remarkably, also in 
this case the calculated VBO is 0.07 eV, although differently divided between AE^ and AV. 
iFmrn our study, however, we cannot fully exclude an effect of ordering on the VBO: we 
can only conclude that this effect, if present, is small and could be hidden by our numerical 
accuracy. 

A further comment is in order concerning the effect on VBO of the possible composition 
fluctuations: we account for them considering that the real composition in the region where 
the band offsets are detected could be sligthly different from the nominal one. In high-quality 
samples this difference should not exceed ~ ±5%, so that the corresponding variation in the 
VBO should be of the order of our numerical accuracy, less than 0.05 eV, estimated from 
the VBO between the end-points. 

Finally, we can estimate the CBO for the heterojunction from the calculated VBO and 
from the experimental data for the gaps of the alloys,— taking into account a quadratic 
interpolation between the end-points: we obtain CBO=— 0.29 eV, with Ino.75Alo.25As higher. 
The CBO is more sensitive than the VBO to the effects of order: although we do not find 
appreciable variations for the VBO, variations up to about 0.14 eV in the CBO could be 
predicted simply by combining the effects of the variations of the band gap in the two alloys 
discussed in the previous Section. 



9 



V. HETEROJUNCTIONS BETWEEN THE (GAAS, ALAS, INAS) BINARIES 



We focus now on the three hmiting cases of the In^^Gai-xAs/In^Ali-yAs hetero junction 
corresponding to interfaces between pure binaries: (a) {x,y) = (0,0), (b) {x,y) = (1,0), (c) 
{x,y) = (0, 1), i.e. GaAs/AlAs, InAs/AlAs and GaAs/InAs interfaces respectively. In these 
simple cases the supercells contain 12 atoms and have a total length equal to approximately 
3ao, depending on the substrate and consequently on the strain state. We are going to 
consider different possible strain configurations within the range of pure GaAs (AlAs) or 
pure InAs substrate. Although the accidental local formation of islands of pure GaAs and 
AlAs binaries at In^Gai_j.As/Inj,Ali_j,As hetero junction is rather unrealistic for the In-rich 
nominal composition that we have considered, it could be eventually possible for InAs. More 
realistically, the cases that we are going to discuss correspond to different heterostructures 
intentionally grown. Since the purpose of this Section is to estimate the possible range of 
variability of the offsets, for the sake of semplicity we will report VBO and CBO values 
without self-energy corrections: this choice does not change our conclusions. 

A. Alloy substrate 

We first consider the case of in-plane lattice constant between those of GaAs (AlAs) and 
InAs. For the sake of definiteness we consider the one of the Ino.75Gao.25As and Ino.75Alo.25As 
alloys, a||=5.86 A. 

Each pure binary semiconductor in this condition is strained: GaAs and AlAs are under 
a tensile in-plane strain (ey =0.058 and 0.047 respectively) and therefore shrink along the 
growth direction, whereas InAs is under a compressive in-plane strain (ey = —0.016) and 
therefore elongates along the growth direction. Considering the theoretical lattice param- 
eters and elastic constants in Table HI we predict a±=5.28, 5.34, 6.07 A and e± = —0.048, 
—0.047, 0.018 for GaAs, AlAs and InAs respectively. A full optimization of the atomic 
positions gives a small overstrain in the interface region. 

From the supercells with optimized atomic positions we calculate the potential lineups, 
on top of which we add the band edges of the binaries in their appropriate strain state 
(see Table HVl and Appendix I A Hi . Finally, the VBO including macroscopic and local strain 
and spin-orbit effects (not self-energy corrections) is: 0.16 eV, 0.30 eV, 0.44 eV for the 
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GaAs/InAs, InAs/AlAs, and GaAs/AlAs strained interfaces (see Fig. panel c)). 

The conduction band bottom edge is also affected by the strain field. Details of the 
calculations are given in Appendix lA 21 distinguishing the cases of direct /indirect gap. Using 
together our numerical calculations and experimental data, the best estimate for the gaps 
of the three binaries strained on the in-plane alloy lattice constant are: E'^g^cfaA^ ~ 0-69 
eV, Egj^'-^s'^ = 0.44 eV and Ef^^^iXs''' = ^-^^ (indirect). /^From the previously calculated 
VBO we can therefore easily predict also the CBO: 0.41 eV, —1.02 eV, and —0.63 eV for 
GaAs/InAs, InAs/AlAs, and GaAs/AlAs hetero junctions respectively. We summarize such 
results again in Fig. El panel c). 

A different kind of band alignment is predicted for the GaAs/AlAs and AlAs/InAs case 
with respect to GaAs/lnAs. Due to the large gap of AlAs, in the former cases the conduction 
and valence band edges of the smaller bandgap material (GaAs or InAs) are both within 
the bandgap of AlAs, that is a type / alignment. At variance, in the latter case, we have a 
type //alignment, with an effective bandgap of the heterostructure of 0.28 eV. 

Finally, we notice that the predicted VBO and CBO for the three interfaces well satisfy 
the transitivity rule within the numerical error, confirming the validity of the linear response 
theory^ for the class of isovalent heterojunctions. 

B. Pure binary substrates 

We have studied also the VBO and the CBO for same three interfaces between the 
binaries, but considering ay = ao(InAs)=5.96 A, (Fig. |3J panel b)), and ay = 5.57 A ^ 
ao(GaAs) ~ ao(AlAs) (Fig. HI panel d)), corresponding respectively to a pure InAs and to a 
pure GaAs or AlAs substrate. 

When a|| = ao(InAs), both GaAs and AlAs are strained. The calculated VBO for 
GaAs/InAs, InAs/AlAs, and GaAs/AlAs is 0.27 eV, 0.07 eV, and 0.34 eV respectively, 
and the corresponding CBO is 0.32 eV, —1.10 eV, and —0.78 eV. 

When a|| = 5.57 A ^ ao(GaAs) ~ ao(AlAs), instead, only InAs shows a sizeable strain. 
The calculated VBO for GaAs/InAs, InAs/AlAs, and GaAs/AlAs is in this case -0.27 eV, 
0.83 eV, and 0.44 eV respectively and the corresponding CBO is 0.71 eV, —0.95 eV, and 
—0.36 eV respectively. We notice that deviations from the transitivity rule up to f^^O.l eV, 
definitely larger than our numerical accuracy, occur for this substrate, and can be ascribed 
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to non neglegible local strain effects at InAs/GaAs(AlAs) interface. 

Summarizing, the main result is that for all the three interfaces the type of alignment 
remains the same by changing ay from ao(InAs) (panel b)) to ao(Ino.75Ga(Al)o.25As) (panel 
c)). At variance, due to sizeale strain effects, when an = 5.57 A ^ ao(GaAs) ~ ao(AlAs) 



(panel d)) the band line-up 



already found in Refs. 



a 



U44 45 



'or the GaAs/InAs interface changes from type //to type I, as 



VI. CONCLUSIONS 



We have estimated VBO=0.11 eV for the In^^Gai-^^As/lnj^Ali-yAs heterojunction with 
homogeneous composition x = y = 0.75, including spin-orbit and self-energy corrections to 
the calclauted LDA value. Neither different possible structural orders nor realistic accidental 
composition fluctuations around the nominal value can modify the VBO more than the 
estimated ^0.1 eV. 

Instead, we have predicted a maximum range of variability of about 1.1 eV for the VBO at 
In3.Gai_a;As/Inj,Ali_j,As heterostructures taking into account interfaces intentionally formed 
between different pure binaries and on different substrates. Neglecting self-energy corrections 
in this case — which could affect a little bit the individual band offset values but are not 
expected to change qualitatively its range of variability — the VBO goes from 0.83 eV in the 
case of InAs/AlAs interface with a parallel lattice constant equal to the one of GaAs or AlAs, 
to 0.44 eV in case of GaAs/AlAs unstrained interface, up to 0.07 eV at InAs/AlAs on InAs 
substrate and —0.27 eV in the case of GaAs/InAs interface on GaAs substrate. The range 
of variation of the CBO is even larger, due to the larger effects of strain on the conduction 
band edges and energy gaps rather then on the valence band edges: the maximum range of 
variation is almost 1.8 eV, from 0.71 eV occurring at GaAs/lnAs interface with a parallel 
lattice constant equal to the one of GaAs, up to —1.10 eV at InAs/AlAs interface with a 
parallel lattice constant equal to the one of InAs. 
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APPENDIX A: STRAIN AND SPIN-ORBIT EFFECTS ON THE BAND EDGES 

1. Valence band edges 

The strain shifts and splits the band edge manifolds of each constituent into the heavy 
hole (HH), light hole (LH) and split-off (SO) states, which can be expressed as: 

EvIhh/LH/SO = K^av + O-vtre + AEy^HH/LH/SO (Al) 

where E^^^ is the average of the valence band top edge manifold at F in the unstrained 
compound, the valence band deformation potential taking into account the effect of the 
hydrostatic component of the strain, and AE^^^^^^^'^ account both for the effects of the 
shear strain and of spin-orbit coupling:— 

AEy^LH = -^Ao + ^<5E,,ooi + ^[Ao' + Ao(5^.,ooi + 7(5£^.,ooi)']^ (A2) 
o 4 z 4 

111 9 1 

AE^ so = — 77A0 + -SEy — - [Aq + AqSE^ ooi + -(5i?„ ooi) ]^ 
o 4 z 4 

where SE^^qqi is the splitting due to the shear strain only, which be obtained directly from 
standard non relativistic band structure calculations of the constituent in the corresponding 
strain state. The final estimate of the valence band levels can be obtained from Eqs. IA2I 
using the calculated SE^^qqi and the experimental spin orbit splitting Aq. 

2. Conduction band edges 

In order to study the effects of strain on the conduction band edges we have to distinguish 
two cases. The conduction band bottom for direct gap semiconductors like GaAs and InAs 
occurs at the T point, it is non degenerate and its position with respect to the reference 
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electrostatic potential depends only to the hydrostatic component of the strain through the 
conduction hand deformation potential';^ 

At variance, in case of AlAs the minimum of the conduction band occurs along the <001> 
direction A, close to the X point, it is six-fold degenerate in the unstrained semiconductor 
since there are six symmetry equivalent <001> directions. Under strain, it is affected both 
by a shift due to the hydrostatic component of the strain through an indirect conduction hand 
deformation potential Oc^indirect and by a splitting due to shear components. In summary, 
and analogously to the expression of Eg. lAll for the valence band top edge manifold, we can 
write for the lowest state of the conduction band bottom: 

EctTtorni^) = E^{A) + acindirecttre + 6Ec ,bottom,001 ■ (A3) 

Similar expressions could be written for the other states of the conduction band manifold. 

In order to give the best estimate for the conduction band offset, due to the limitation 
of the LDA, we can use the equations above without spin-orbit to extract from the non- 
relativistic calculations of the materials in their strained state the quantities SE^^ooi and 
SEc^hottomfloi, and then we can use them again by inserting the experimental values of all the 
other band parameters {Eg, ay, ac, Ocmdirect) for each material to get: 

^j^strained ^ ^0 ^ _ ^^^^^^ _ AEy,HH/LH (A4) 

where LH and HH holds respectively for GaAs and for InAs (their topmost valence states, 
see Table HV|) . and 

Eg^Jndirect = ^g.mdirerf + {'^c,indirect " ay)tre + 6 E ^fiOl , bottom — AE^^ffH (A5) 

for AlAs. Here Eg = E*^ — E^ + ^ (experimental unstrained gap, without spin-orbit), and 
5E,,ooi,bottom < and AE^^/^^ > 0. 

APPENDIX B: ROLE OF d ELECTRONS 

It has been shown^ that the d electrons affect the properties of the II- VI zinc-blende 
compounds through p — ci repulsion and hybridization. In particular, they lead to a reduction 
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of spin orbit splitting and increase of the valence band offset between common anion systems 
(by pushing up the Valence Band Maximum) with respect to calculations where d electrons 
are considered as frozen (core) states. 

In 111-V systems, the d electrons should play a minor role. In this Appendix we discuss in 
particular this point for our systems of interest by performing additional first-principle all- 
electron calculations on the bulk binary systems and their interfaces both considering Ga(In) 
d electrons as core states or as relaxed (valence) states and we compare these results with 
the pseudopotential calculations. To this aim we use a full-potential linearized augmented- 
plane-wave code (FLAPW)^2i^ using DFT-LDA, where spin-orbit splitting can be included 
as a perturbation on the semirelativistic calculation. This allows us to discuss the effect 
of d electrons also on spin-orbit. We will consider the effects on the band structure of 
free-standing and strained GaAs and InAs bulks. 

The calculated FLAPW equilibrium lattice constant {a.eq) for GaAs and InAs are 
5.68(5.63) A and 6.07(6.02) A respectively when d electrons are considered as core (va- 
lence) states. The slight increase of the lattice constant when d electrons are frozen is along 
the trend of similar calculations . 

The calculated spin-orbit splitting (Aq) does not change if d electrons are treated as core 
or valence states, providing that the corresponding equilibrium theoretical lattice parameter 
is considered: it is equal to 0.36 eV for GaAs and 0.39 eV for InAs (the same as the 
experimental one). Therefore we can conclude that the influence of d states on the spin- 
orbit splitting of valence band maximum is negligible. This justifies our choice of using a 
posteriori the experimental values of the spin-orbit splittings to correct the non-relativistic 
LDA pseudopotential structure. 

We have also studied the role of d states when strain is applied, considering the same strain 
states used in our pseudopotential calculations both for GaAs strained over InAs (ay =5.96 
A, see Table HV|l and for InAs strained over Gao.5Alo.5As (a||=5.57 A, see Table HV|l . We 
choose to keep the same strain tensor in FLAPW and PWSCF calculations for a systematic 
comparison of the results, although the two different computational method would give 
slightly different equilibrium lattice constants for the binary compounds and consequently 
slightly different strain states. In these cases, we have used the elastic constants reported in 
Table lU a slightly different value of as calculated with the FLAPW method should lead 
to a negligible variation of the ratio The calculated 5E^^q{ is equal to —0.17 (—0.14) 
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eV and SE^^l'^l= -0.35 (-0.35 ) eV for d states treated as core (valence). Therefore, a small 
variation in the strain splitting effect is possible due to a different treatment of d states, but 
very limited, here estimated within 0.03 eV. 
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TABLE I: Relevant structural and electronic equilibrium parameters for bulk GaAs, AlAs, InAs. 
For AlAs we report the indirect /direct energy gap. Experimental data are reported in round 
bracketsii 





GaAs 


AlAs 


InAs 


ao (A) 


5.55 (5.65) 


5.60 (5.66) 


5.96 (6.06) 


B{kBar) 


760 (784) 


750 (773) 


670 (579) 


cii{kBar) 


1240 (1221) 


1130 (1250) 


940 (833) 


ci2{kBar) 


520 (566) 


560 (534) 


540 (453) 


EgieV) 


1.44 (1.52) 


1.35 (2.24) / 2.41 (3.10) 


0.39 (0.42) 
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TABLE II: Calculated and experimental mismatches between binary semiconductors. 



Mismatch (%) 


Theory 


Experiment 


GaAs-AlAs 


1.0 


0.1 


GaAs-InAs 


7.5 


7.2 


AlAs-InAs 


6.4 


7.1 
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TABLE III: Nearest-neighbour (NN) and next-nearest-neighbour distances (NNN) in 
Ino.75Gao.25 As and in Ino.75Alo.25As alloys in the three ordered structure considered in the text. 
We indicate with cat.-cat. all the possible pairs of cations. The error in the distances indicates 
their spreading in the structure considered. Units are A. 





S 


L 


F 


Ino.75Gao.25 As 


NN(In-As) 
NN(Ga-As) 


2.57±0.01 
2.45 


2.58 
2.43 


2.59±0.01 
2.41 


NNN {cat-cat.) 
NNN (As- As) 


4.15±0.07 
4.07±0.16 


4.15 
4.15±0.18 


4.15 
4.146±0.22 


Ino.75Alo.25As 


NN(In-As) 

NN(Al-As) 


2.58±0.01 
2.45 


2.58 

2.44 


2.58±0.01 

2.43 


NNNicat.-cat.) 
NNN(An-An) 


4.15±0.07 
4.083±0.15 


4.15 
4.15±0.16 


4.15 
4.15±0.19 
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TABLE IV: Calculated structural and band parameters for the three binaries GaAs, AlAs, InAs in 
three pseudomorphic structures: with a parallel lattice constant ay =5.96 A (equal to that of bulk 
InAs), a||=5.86 A (bulk Ino.75Ga(Al)o.25As), a||=5.57 A (bulk GaAs or AlAs). For the energy gap, 
the values reported in this table arc those for the minimum gap only (indirect in AlAs), and are 
obtained by taking from self-consistent calculations the effects of strain and by adding a posteriori 
an empirical correction from the comparison of experimental/theoretical data for the binary bulks. 





GaAs 


AlAs 


InAs 


a||=5.96 A 






0.075 


0.064 


0.000 






-0.064 


-0.063 


0.000 


a± (A) 




5.19 


5.25 


5.96 


^E^,LH/HH/SO 


(eV) 


0.44 (LH) 


0.36 (LH) 


0.13 (HH,LH) 


E, (eV) 




0.47 


1.59 


0.42 


a||=5.86 A 






0.058 


0.047 


-0.016 


e± 




-0.048 


-0.047 


0.018 


a± (A) 




5.28 


5.34 


6.07 




(eV) 


0.37 (LH) 


0.29 (LH) 


0.19 (HH) 


E, (eV) 




0.69 


1.76 


0.44 






a =5.57 


A 




ey 




0.005 


-0.005 


-0.065 


e± 




-0.004 


0.005 


0.075 


a± (A) 




5.52 


5.63 


6.41 


^E^,LH/HH/SO 


(eV) 


0.13(LH) 


0.11 (HH) 


0.41 (HH) 


E, (eV) 




1.45 


2.25 


0.47 
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FIG. 1: Band structure of Ino.75Gao.25As (top panels) and Ino.75Alo.25As (bottom panels) alloys in 
the (001)3+1 superlattice (S), luzonite (L) and famatinite (F) structures. Thick lines emphasize 
valence and conduction band edges. The zero of the energy scale is always set at the valence band 
maximum. 

FIG. 2: Density of states of Ino.75Gao.25As and Ino.75Alo.25As alloys in the (001)3+1 superlattice 
(S), luzonite (L) and famatinite (F) structures. The zero of the energy scale is always set at the 
valence band maximum. 

FIG. 3: Macroscopic average^ of the charge density pMacro (solid lines, right scale) and electrostatic 
potential VMacro (dashed lines, left scale) for the Ino.75Gao.25As/Ino.75Alo.25As heteroj unction when 
the alloys are described using a luzonite structure (upper panel) and a (001)3+1 superlattice (lower 
panel). The position of the anionic planes is indicated in the x-axis. 

FIG. 4: Calculated valence (conduction) band offsets for the (001) interfaces between the alloys 
(top panel a)) and binaries GaAs, AlAs and InAs pseudomorphically grown on a substrate with: 
b) a||=5.96 A (equal to that of bulk InAs), c) a||=5.86 A (bulk Ino.75Ga(Al)o.25As), d) a||=5.57 A 
(bulk GaAs or AlAs). The band alignments are calculated between and Ec, the highest valence 
(lowest conduction) states including spin-orbit effects. We indicate in boldface the unstrained 
compounds. Values reported here include spin-orbit but not self-energy corrections. 
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